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We derive the parallel upper critical field, H C 2, as a function of the temperature T in quasi-2D 
organic compound A-(BETS)2FeCl4, accounting for the formation of the nonuniform LOFF state. 
To further check the 2D LOFF model we propose to study the H C 2(T) curve at low T in tilted fields, 
where the vortex state is described by the high Landau level functions characterized by the index 
n. We predict a cascade of first order transitions between vortex phases with different n, between 
phases with different types of the symmetry at given n and the change of the superconducting 
transition from the second order to the first order as FeCU ions are replaced partly by GaCk ions. 
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Very recently magnetic-field-induced superconductiv- 
ity has been observed in the quasi two-dimensional (2D) 
organic conductor A-(BETS)2 FeCU @- At zero mag- 
netic field the antiferromagnetic ordering of Fe 3+ mo- 
ments in this compound gives rise to a metal-insulator 
transition at temperatures around 8 K. A magnetic field 
above 10 T restores the metal phase. Further increase 
of the magnetic field induces the superconducting tran- 
sition at B ss 17 T. As it has been revealed in the very 
high magnetic field experiments 0, the maximum criti- 
cal temperature T c w 4 K is reached at Bg ss 33 T and 
T C (B) drops to 2 K at B = 42 T. Such an unusual be- 
havior was interpreted in Refs. jj], 0] as a manifestation 
of the Jaccarino-Peter (JP) effect M, when the exchange 
field of aligned Fe 3+ spins compensates the external field 
in their combined effect on the electron spins. The JP 
compensation effect has been already proved to be re- 
sponsible for the magnetic-field-induced superconductiv- 
ity in pseudo-ternary Eu-Sn molybdenum chalcogenides 
J§ and in CePb 3 §]. In contrast to A-(BETS) 2 FeCl 4 , 
in Eu-Sn molybdenum chalcogenides the superconduct- 
ing transition also exists in the absence of a magnetic 
field, and as the magnetic field increases, it first destroys 
superconductivity, but further restores it in the range 
of fields 4 T< B < 22 T at very low temperatures j|. 
In CePb3 in zero magnetic field the heavy-fermion anti- 
ferromagnetic ground state was found g|. Similarly to 
A-(BETS)2FeCi4, the magnetic field destroys the anti- 
ferromagnetic ordering and induces superconductivity in 
fields higher than 14 T. 

High-field superconductivity exists in all these com- 
pounds not only due to the JP effect, but also because 
of strong suppression of the orbital effect of the external 
field. In CePb 3 and Eu-Sn molybdenum chalcogenides 
the orbital effect is weak due to a very short supercon- 
ducting coherence length in the presence of impurities. 
In contrast, in the clean A-(BETS)2FeCLt superconduc- 
tor strong quasi-2D anisotropy leads to negligible orbital 
effect for the field applied parallel to the layers. Here the 
parallel upper critical field is controlled by the spin ef- 
fect only. Then the temperature dependence of H C 2 may 



be well described by the quasi-2D (Josephson-coupled- 
layers) BCS model [||. In this model at low tempera- 
tures T < 0.55 T c , the Larkin-Ovchinnikov-Fulde-Ferrell 
(LOFF) state is a correct superconducting phase |). 
We note that previous attempts to identify unambigu- 
ously the LOFF state were unsuccessful because in 3D 
crystals orbital effect dominates over spin effect, while 
quasi-2D systems, like intercalated dichalcogenides, were 
not clean enough to allow for the LOFF state. The very 
clean 2D organic crystal A-(BETS)2FeCi4 with clear evi- 
dence for dominant spin effect is the most favorable sys- 
tem to detect the LOFF state unambiguously. 

In this Letter we use a 2D LOFF model for the clean 
superconductor A-(BETS) 2 FeCl 4 to calculate the H c2 (T) 
curve for the orientation of the magnetic field parallel to 
the layers, accounting for the exchange field induced by 
Fe 3+ polarized magnetic moments. We show that ac- 
counting for the 2D LOFF state improves fitting at low 
temperatures in comparison with the 3D model as was 
also discussed in Ref. In addition to results 0, we 
calculate H C 2 for the magnetic fields tilted by angle 9 
with respect to the layers. We show that a quasiperiodic- 
like angular dependence of H C 2 is a specific property of 
the vortex state, which originates from the 2D LOFF 
phase when the orbital effect, caused by the perpendic- 
ular component of the external field, is turned on as 6 
increases. This vortex structure is characterized by the 
discrete variable n, the Landau level index. Changes of 
n with 8 lead to a cascade of first order transitions of the 
vortex lattice. We argue also that partial replacement 
of FeCU ions by nonmagnetic GaCU ions should lead to 
the suppression of the LOFF state and to the first order 
normal-to-superconducting (N-S) transition. 

We consider first the H C 2 (T) curve for the parallel ori- 
entation of the external field. The total effective field 
acting on the electron spins is H e f f = B — I (S) / fi, where 
fj, « [lb is the electron magnetic moment, and I is the 
exchange integral for the interaction between the conduc- 
tion electrons and Fe 3+ magnetic moments. The negative 
sign implies the JP compensation effect. In the whole 
(B, T)-region, where the field-induced superconductivity 
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exists in A-(BETS) 2 FeCl 4 , the Fe 3+ moments are satu- 
rated, (S) = Sq. The maximum T C (H) corresponds to 
H e ff = 0. Hence, the effective field acting on spins is 
simply H e ff = B — Bq, where Bq = ISq/h. The up- 
per critical field is given by the solution of the linearized 
equation for the superconducting order parameter A(r), 
where r = (x, y) is the in-plane coordinate: 



T„=4.3 K 



A(r) = J iv" (r-r')A(r')rfr', 



(1) 



For the nonuniform order parameter, A(r) = 
Aexp(iq • r), the kernel Kq(t — r') takes the form || |): 



Jfo(q) = 27r|A|TRe^ I 



[27r|v F (Z) |] — 
lo + ifj,H e ff + iwp(t)<\/2 ' 



(2) 



where u> denotes the Matsubara's frequencies, A is the 
coupling parameter, line integral over the coordinate I 
is along the 2D Fermi surface in the momentum space, 
Vf(1) is the in-plane Fermi velocity. 

As has been demonstrated by Larkin and Ovchinnikov 
and Fulde and Ferrell 0, in magnetic fields above some 
critical value the modulated superconducting phase with 
nonzero q gives higher T c for singlet pairing. A uniform 
superconducting order parameter in the presence of uni- 
form Zeeman splitting is unfavorable because electrons 
with opposite spins and momenta have different ener- 
gies e(fc) ± fj,H, where e(k) is the electron kinetic energy. 
Meanwhile, the Cooper instability towards formation of 
superconducting pairs is strongest when the energies of 
electrons in a pair are close to each other. Hence, pairing 
of electrons with momenta k ± q at q ~ fiH e tf/vp and 
opposite spins is more favorable than at q = 0, because 
energies e(k + q) — tp — \iH e f f and e(k — q) — ep + \iB e ff 
are closer. Obviously, gain in the superconducting energy 
depends on the dimensionality of the electron motion. In 
a quasi ID system one can compensate Zeeman splitting 
almost completely, while in the 2D case compensation 
is only partial due to different orientations of the Fermi 
velocity, and in 3D crystals compensation is even less ef- 
fective. Hence, H C 2 is very sensitive to the geometry of 
the in-plane Fermi surface (FS). 

In the following we consider an isotropic in-plane elec- 
tron spectrum and also the anisotropic one defined as 



e(p) = W{cosp x + Tcospy), e F = 0, 



(3) 



where the parameter T < 1 accounts for the in-plane 
anisotropy. T C (B) for such a nonuniform state is defined 
by the equation K(q,B,T) = 1. The wave- vector q at 
the transition from the normal state to the LOFF state is 
that which gives a maximal H C 2 at a given T. It depends 
on the vp value, but H C 2 does not. The corresponding 
H C 2(T) = B ± H e ff curves for the 3D spherical FS, 
2D circular FS and for the spectrum given by Eq. (|3|) at 
T = 0.7, calculated for T c = 4.3 K, are shown in Fig. |l| 
together with the experimental data 0. 
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FIG. 1: The dependence of the critical temperature on the 
magnetic field parallel to the layers in the presence of the 
Jaccarino-Peter effect. The thin solid curve is the transition 
from the normal metal to the uniform superconducting state. 
The dotted (thick) line is the transition to the LOFF state in 
the isotropic 3D (2D) system and the dashed line is for the 
in-plane anisotropic 2D system. Dots are experimental data 
obtained by resistivity measurements [2] . 



The common part of the B C 2 (T) curves corresponding 
to the transition into the uniform state is presented by 
the thin solid line. To give insight into the importance 
of the in-plane anisotropy we calculate the critical value 
of h = 2{iH e f f I Ao at T = as a function of T. Here 
Ao = 1 .76T C is the superconducting gap at T = 0. The 
upper critical value of h is given by the equation 



^/ 2 lnl^-g^i-^cos 2 !*)! 
du — ■ ■ 

i) v i — r 2 cos 2 u 



= 0. 



(4) 



At r = 0.7 we obtain the upper critical value h cr = 2.64 
corresponding to ^H e ff = f.32A . For comparison, we 
obtain h cr — 2, corresponding to [iB e ff = Ao, in the 
case of a circular FS. In this approach h cr (T) increases 
without saturation as T decreases. 

To check whether the LOFF state is indeed realized in 
A-(BETS)2FeCi4 a decisive experimental test is needed. 
Here we propose that such a critical experiment may be 
the study of the dependence H C 2(T) for tilted-toward- 
the-plane field orientation. The perpendicular compo- 
nent of the magnetic field suppresses pairing by the or- 
bital mechanism and leads to the formation of the vor- 
tex state on the already nonuniform LOFF background 
caused by the spin effect. The interplay between orbital 
and the spin effects gives rise to a very peculiar upper- 
critical-field behavior || ||] resulting from the solutions 
with higher Landau level functions for the superconduct- 
ing order parameter in the vortex state. 

To derive the H C 2(T) curve at tilting angles 9 we ac- 
count for the effect of the vector potential describing the 
perpendicular component of the magnetic field, B^, in 
addition to the spin effect of -By and of the exchange 
field. Now, in Eq. (1) the kernel is given as 



K(r, r') = Kq(y — r') exp 



2ie/hc / A(s)ds 



(5) 
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FIG. 2: The dependence of the upper critical magnetic field 
on T at different tilt angles 6 for an isotropic 2D Fermi surface. 
Different parts of curve correspond to different Landau level 
functions characterized by the index n. 



for an isotropic 2D FS and the in-plane vector-potential 
is A = — (Bj_ x r)/2, with B ± = 5sin6», @. Following 
Refs. H |j| general solutions of the integral equation for 
the order parameter are the Landau level functions: 



A„(r) = exp(-in</? - p 2 /2)p v 



(6) 



where p = r^J eB±/hc and (p is the polar angle. 

When the orbital effect dominates, it is the zero Lan- 
dau level function which gives the solution corresponding 
to the highest magnetic critical field. However, if the spin 
effect is strong and the electron mean free path is much 
bigger than the superconducting correlation length (i.e. 
in the clean limit), the higher Landau level functions for 
a superconducting nucleation give higher H C 2, because 
they provide rapid spatial variations of the supercon- 
ductor order parameter needed to compensate Zeeman 
splitting of electrons in a singlet pair. Discontinuous ad- 
justments of the discrete variable n to the change of the 
angle 9 give rise to unusual quasi-oscillatory angular and 
temperature dependence of H c2 , see Refs. ]|, [|. The 
general expression for ff c2 is given by the equation 



ln(T/T c ) = ttT x 
Re H {/ 



(7) 



(-l) n L n (x) exp(-a;/2)dx _ 2 
[(lo + ipH eff ) 2 + xv 2 F eB±/4:hc} V 2 uj 



where L n (x) are the Laguerre polynomials. In Fig. [2] we 
show the expected behavior of H C 2 for different tilting 
angles 9 for the model with isotropic 2D FS. In these 
calculations we take the value 5 T for the orbital magnetic 
critical field at T = as was obtained for (BETS^GaCU 
with similar conducting layers jlj. The behavior of the 
system in the quasi-lD limit may be calculated in the 
same way basing on the results of Ref . . 

Next we discuss the vortex lattice originating from the 
LOFF state due to the orbital effect of B^. Inside the 
superconducting phase, the order parameter must form 
a lattice and without the spin effect it is a triangular 
Abrikosov vortex lattice made of n — Landau level 
functions. In such a standard vortex lattice the number 



of zeros of the superconducting order parameter coincides 
with number of quantized vortices, i.e. there is one zero 
per unit cell. However, for higher Landau level solutions 
the number of zeros exceeds number of quantized vortices 
determined by the component B± and only a fraction 
of the zeros coincides with the centers of quantized vor- 
tices. Respectively, zeros become inequivalent, they form 
a lattice determined by two length scales, the average 
distance between quantized vortices, a = (Qq/B^) 1 / 2 , 
and the LOFF period determined by the effective field, 
I = hvp/pH e ff. This leads to a large variety of exotic 
vortex lattices as has been predicted in Refs. |u| [l4|. 

To determine the equilibrium vortex structure we 
minimize the free energy. In Ref. Jl3| the general- 
ized Ginzburg-Landau (GL) functional has been derived 
to describe the LOFF state near the tricritical point, 
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1.07 T c /fi, T* = 0.55 T c , for isotropic FS: 
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0.85 £ 2 |A| 2 |VAf 



(A*VA) 2 + (AVA*) 2 



where N(0) is the DOS at the Fermi level, V = V - 
(2ie/hc)A and £0 = hup I2txT c is the superconducting 
correlation length. Note that we need to include in the 
functional the second-order derivative term ~ |V 2 A| 2 
due to the negative sign of the first-derivative term. Sim- 
ilarly, we add the sixth order term ~ |A| 6 Strictly 
speaking the derived functional is valid near the tricriti- 
cal point only, but we may expect that it describes qual- 
itatively correctly (as the standard GL functional) the 
LOFF state at all temperatures. Extending the analysis 
of the vortex lattice structure jl2| to this new functional, 
one can calculate the structures of possible new vortex 
lattices. It turns out that the N-S transition becomes of 
the first order in some temperature and angle intervals 
|l3| and the lattice structure exhibits zeros with posi- 
tive and negative winding numbers w, \w\ > 1, as well 
as strongly anisotropic structures. In Fig. [| we present 
as an example the square asymmetric vortex lattice ob- 
tained numerically by use of the functional (]8|). It cor- 
responds to the Landau level function n = 2 and have 
five zeros in the unit cell, two of them has w = —1, while 
three have w = +1, so that there is one flux quantum per 
unit cell. Note that zeros are positioned inequivalently 
and that different types of symmetry in their positions 
as well as different indecies n lead to a variety of vortex 
lattices and first order phase transitions between them. 

Let us discuss now the anticipated behavior of the 
system A — (BETS)2(FeCl4)i_ K (GaCl4)a; using informa- 
tion |l(| on H c2 (T) at x — 1. Replacing the Fe 3+ ion 
randomly by nonmagnetic Ga 3+ reduces the average ex- 
change field as Bq(1 — x). Additionally, a random ex- 
change field leads to electron scattering with rate of or- 
der pBox. It causes suppression of the LOFF state at 
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FIG. 3: Square asymmetric vortex lattice described by the 
Landau level functions n = 2 formed in the presence of or- 
bital and spin effects. The profiles of the amplitude of the 
superconducting order parameter are shown. The dashed line 
represents the unit cell with five zeros. Their winding num- 
bers, ±1, are shown. There is one flux quantum per unit 
cell. 

x ~ 1/2, when the scattering rate becomes comparable 
with Ao- Hence, as x increases, while remaining small, 
in Figs. |l|, ||the curves T(H) scale as (1 — x)H, while the 
thick and dashed lines approach the dotted line due to the 
suppression of the LOFF state. At x w 1/2 we anticipate 
that due to strong exchange scattering the LOFF state 
vanishes, and the first order superconducting transition 



takes place near the dotted line, at H e ff — Aq/^V^ at 
low T. Observation of the change of the transition type 
with x would confirm unambiguously existence of the 
LOFF state at x = 0. Interestingly, as x tends to unity, 
the system becomes clean again, and the LOFF state may 
reappear. In fact, at x = 1 the value of H C 2 ~ 14 T at 
low T is close to that anticipated for the 2D LOFF state. 
We note that T c increases slightly from 4.3 K at x = to 
5.5 K at x = 1. Additionally, anisotropy probably drops 
with x, but even at x = 1 the system remains close to 
2D, because the inter layer correlation length is about the 
same as the interlayer spacing [[t6|. Hence, both effects 
do not change our qualitative predictions. 

In conclusion, in A — (BETS^FeCLj we anticipate 
first order phase transitions between the superconduct- 
ing phases corresponding to different higher Landau level 
states and between lattices with different symmetry. We 
anticipate also the change of the S-N transition from sec- 
ond order to first order as x approaches 1/2 in the system 
A - (BETS) 2 (FeCl4)i-x(GaCl 4 ) 2: . These phase transi- 
tions may be detected by magnetization and specific heat 
measurements. We think that neutron scattering, fiSK, 
and NMR measurements may be used additionally to re- 
veal the structure of peculiar vortex phases formed due 
to the interplay of competing spin and orbital effects in- 
duced by the external and exchange fields. 
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